Aedes albopictus colonies from different geographic origins differ in their sleep and activity levels but not in the time of peak activity.

Mosquitoes occupy a wide range of habitats where they experience various environmental conditions. The ability of some species, such as the tiger mosquito, Aedes albopictus, to adapt to local conditions certainly contributes to their invasive success. Among traits that remain to be examined, mosquitoes' ability to time their activity with that of the local host population has been suggested to be of significant epidemiological importance. However, whether different populations display heritable differences in their chronotype has not been examined. Here, we compared laboratory strains originating from 8 populations from 3 continents, monitored their spontaneous locomotor activity patterns, and analyzed their sleep-like states. Overall, all strains showed conserved diurnal activity concentrated in the hours preceding the crepuscule. Similarly, they all showed increased sleep levels during the morning and night hours. However, we observed strain-specific differences in the activity levels at each phase of the day. We also observed differences in the fraction of time that each strain spends in a sleep-like state, explained by variations in the sleep architecture across strains. Human population density and the latitude of the site of geographic origin of the tested strain showed significant effects on sleep and activity patterns. Altogether, these results suggest that Ae. albopictus mosquitoes adapt to local environmental conditions via heritable adaptations of their chronotype.

A genotyping array for the globally invasive vector mosquito, Aedes albopictus.

BACKGROUND: Although whole-genome sequencing (WGS) is the preferred genotyping method for most genomic analyses, limitations are often experienced when studying genomes characterized by a high percentage of repetitive elements, high linkage, and recombination deserts. The Asian tiger mosquito (Aedes albopictus), for example, has a genome comprising up to 72% repetitive elements, and therefore we set out to develop a single-nucleotide polymorphism (SNP) chip to be more cost-effective. Aedes albopictus is an invasive species originating from Southeast Asia that has recently spread around the world and is a vector for many human diseases. Developing an accessible genotyping platform is essential in advancing biological control methods and understanding the population dynamics of this pest species, with significant implications for public health. METHODS: We designed a SNP chip for Ae. albopictus (Aealbo chip) based on approximately 2.7 million SNPs identified using WGS data from 819 worldwide samples. We validated the chip using laboratory single-pair crosses, comparing technical replicates, and comparing genotypes of samples genotyped by WGS and the SNP chip. We then used the chip for a population genomic analysis of 237 samples from 28 sites in the native range to evaluate its usefulness in describing patterns of genomic variation and tracing the origins of invasions. RESULTS: Probes on the Aealbo chip targeted 175,396 SNPs in coding and non-coding regions across all three chromosomes, with a density of 102 SNPs per 1 Mb window, and at least one SNP in each of the 17,461 protein-coding genes. Overall, 70% of the probes captured the genetic variation. Segregation analysis found that 98% of the SNPs followed expectations of single-copy Mendelian genes. Comparisons with WGS indicated that sites with genotype disagreements were mostly heterozygotes at loci with WGS read depth < 20, while there was near complete agreement with WGS read depths > 20, indicating that the chip more accurately detects heterozygotes than low-coverage WGS. Sample sizes did not affect the accuracy of the SNP chip genotype calls. Ancestry analyses identified four to five genetic clusters in the native range with various levels of admixture. CONCLUSIONS: The Aealbo chip is highly accurate, is concordant with genotypes from WGS with high sequence coverage, and may be more accurate than low-coverage WGS.

Population-specific responses to developmental temperature in the arboviral vector Aedes albopictus: Implications for climate change.

The increase of environmental temperature due to current global warming is not only favouring the expansion of the distribution range of many insect species, but it is also changing their phenology. Insect phenology is tightly linked to developmental timing, which is regulated by environmental temperatures. However, the degree to which the effects of developmental temperatures extend across developmental stages and their inter-stage relationships have not been thoroughly quantified in mosquitoes. Here, we used the mosquito Aedes albopictus, which is an aggressive invasive species and an arboviral vector, to study how developmental temperature influences fitness across developmental stages, thermal traits, energy reserves, transcriptome and Wolbachia prevalence in laboratory-reared populations originally collected from either temperate or tropical regions. We show that hatchability, larval and pupal viability and developmental speed are strongly influenced by temperature, and these effects extend to wing length, body mass, longevity and content of water, protein and lipids in adults in a population-specific manner. On the contrary, neither adult thermal preference nor heat resistance significantly change with temperature. Wolbachia density was generally lower in adult mosquitoes reared at 18°C than at other tested temperatures, and transcriptome analysis showed enrichment for functions linked to stress responses (i.e. cuticle proteins and chitin, cytochrome p450 and heat shock proteins) in mosquitoes reared at both 18 and 32°C. Our data showed an overall reduced vector fitness performance when mosquitoes were reared at 32°C, and the absence of isomorphy in the relationship between developmental stages and temperature in the laboratory population deriving from larvae collected in northern Italy. Altogether, these results have important implications for reliable model projections of the invasion potentials of Ae. albopictus and its epidemiological impact.

Virome and nrEVEome diversity of Aedes albopictus mosquitoes from La Reunion Island and China.

BACKGROUND: Aedes albopictus is a public health threat for its worldwide spread and ability to transmit arboviruses. Understanding mechanisms of mosquito immunity can provide new tools to control arbovirus spread. The genomes of Aedes mosquitoes contain hundreds of nonretroviral endogenous viral elements (nrEVEs), which are enriched in piRNA clusters and produce piRNAs, with the potential to target cognate viruses. Recently, one nrEVE was shown to limit cognate viral infection through nrEVE-derived piRNAs. These findings suggest that nrEVEs constitute an archive of past viral infection and that the landscape of viral integrations may be variable across populations depending on their viral exposure. METHODS: We used bioinformatics and molecular approaches to identify known and novel (i.e. absent in the reference genome) viral integrations in the genome of wild collected Aedes albopictus mosquitoes and characterize their virome. RESULTS: We showed that the landscape of viral integrations is dynamic with seven novel viral integrations being characterized, but does not correlate with the virome, which includes both viral species known and unknown to infect mosquitoes. However, the small RNA coverage profile of nrEVEs and the viral genomic contigs we identified confirmed an interaction among these elements and the piRNA and siRNA pathways in mosquitoes. CONCLUSIONS: Mosquitoes nrEVEs have been recently described as a new form of heritable, sequence-specific mechanism of antiviral immunity. Our results contribute to understanding the dynamic distribution of nrEVEs in the genomes of wild Ae. albopictus and their interaction with mosquito viruses.

Computational Methods for the Discovery and Annotation of Viral Integrations.

The transfer of genetic material between viruses and eukaryotic cells is pervasive. Somatic integrations of DNA viruses and retroviruses have been linked to persistent viral infection and genotoxic effects. Integrations into germline cells, referred to as Endogenous Viral Elements (EVEs), can be co-opted for host functions. Besides DNA viruses and retroviruses, EVEs can also derive from nonretroviral RNA viruses, which have often been observed in piRNA clusters. Here, we describe a bioinformatic framework to annotate EVEs in a genome assembly, study their widespread occurrence and polymorphism and identify sample-specific viral integrations using whole genome sequencing data.

Thermal biology of invasive Aedes mosquitoes in the context of climate change.

The increasing incidence of arboviral diseases in tropical endemic areas and their emergence in new temperate countries is one of the most important challenges that Public Health agencies are currently facing. Because mosquitoes are poikilotherms, shifts in temperature influence physiological functions besides egg viability. These traits impact not only vector density, but also their interaction with their hosts and arboviruses. As such the relationship among mosquitoes, arboviral diseases and temperature is complex. Here, we summarize current knowledge on the thermal biology of Aedes invasive mosquitoes, highlighting differences among species. We also emphasize the need to expand knowledge on the variability in thermal sensitivity across populations within a species, especially in light of climate change that encompasses increase not only in mean environmental temperature but also in the frequency of hot and cold snaps. Finally, we suggest a novel experimental approach to investigate the molecular architecture of thermal adaptation in mosquitoes.

Endogenous viral elements in mosquito genomes: current knowledge and outstanding questions.

Integrations from non-retroviral RNA viruses (nrEVEs) have been identified across several taxa, including mosquitoes. Amongst all Culicinae species, the viral vectors Aedes aegypti and Aedes albopictus stand out for their high number of nrEVEs. In addition, Aedes nrEVEs are enriched in piRNA clusters and generate piRNAs that can silence incoming viral genomes. As such, nrEVEs represent a new form of inherited antiviral immunity. To propel this discovery into novel transmission-blocking vector control strategies, a deeper understanding of nrEVE biology and evolution is essential because differences in the landscape of nrEVEs have been identified in wild-caught mosquitoes, the piRNA profile of nrEVEs is not homogeneous and nrEVEs outside piRNA clusters exist and are expressed at the mRNA level. Here we summarise current knowledge on nrEVEs in mosquitoes and we point out the many unanswered questions and potentials of these genomic elements.

Cross Talk between Viruses and Insect Cells Cytoskeleton.

Viruses are excellent manipulators of host cellular machinery, behavior, and life cycle, with the host cell cytoskeleton being a primordial viral target. Viruses infecting insects generally enter host cells through clathrin-mediated endocytosis or membrane fusion mechanisms followed by transport of the viral particles to the corresponding replication sites. After viral replication, the viral progeny egresses toward adjacent cells and reaches the different target tissues. Throughout all these steps, actin and tubulin re-arrangements are driven by viruses. The mechanisms used by viruses to manipulate the insect host cytoskeleton are well documented in the case of alphabaculoviruses infecting Lepidoptera hosts and plant viruses infecting Hemiptera vectors, but they are not well studied in case of other insect-virus systems such as arboviruses-mosquito vectors. Here, we summarize the available knowledge on how viruses manipulate the insect host cell cytoskeleton, and we emphasize the primordial role of cytoskeleton components in insect virus motility and the need to expand the study of this interaction.

Profile of Small RNAs, vDNA Forms and Viral Integrations in Late Chikungunya Virus Infection of Aedes albopictus Mosquitoes.

The Asian tiger mosquito Aedes albopictus is contributing to the (re)-emergence of Chikungunya virus (CHIKV). To gain insights into the molecular underpinning of viral persistence, which renders a mosquito a life-long vector, we coupled small RNA and whole genome sequencing approaches on carcasses and ovaries of mosquitoes sampled 14 days post CHIKV infection and investigated the profile of small RNAs and the presence of vDNA fragments. Since Aedes genomes harbor nonretroviral Endogenous Viral Elements (nrEVEs) which confers tolerance to cognate viral infections in ovaries, we also tested whether nrEVEs are formed after CHIKV infection. We show that while small interfering (si)RNAs are evenly distributed along the full viral genome, PIWI-interacting (pi)RNAs mostly arise from a ~1000 bp window, from which a unique vDNA fragment is identified. CHIKV infection does not result in the formation of new nrEVEs, but piRNAs derived from existing nrEVEs correlate with differential expression of an endogenous transcript. These results demonstrate that all three RNAi pathways contribute to the homeostasis during the late stage of CHIKV infection, but in different ways, ranging from directly targeting the viral sequence to regulating the expression of mosquito transcripts and expand the role of nrEVEs beyond immunity against cognate viruses.

Exploring Changes in the Microbiota of Aedes albopictus: Comparison Among Breeding Site Water, Larvae, and Adults.

The mosquito body hosts highly diverse microbes, which influence different physiological traits of both larvae and adults. The composition of adult mosquito microbiota is tightly linked to that of larvae, which are aquatic and feed on organic detritus, algae and prokaryotic microorganisms present in their breeding sites. Unraveling the ecological features of larval habitats that shape the structure of bacterial communities and their interactions with the mosquito host is still a poorly investigated topic in the Asian tiger mosquito Aedes albopictus, a highly invasive species that is vector of numerous arboviruses, including Dengue, Chikungunya, and Zika viruses. In this study, we investigated the composition of the bacterial community present in the water from a natural larval breeding site in which we separately reared wild-collected larvae and hatched eggs of the Foshan reference laboratory strain. Using sequence analysis of bacterial 16S rRNA gene amplicons, we comparatively analyzed the microbiota of the larvae and that of adult mosquitoes, deriving information about the relative impact of the breeding site water on shaping mosquito microbiota. We observed a higher bacterial diversity in breeding site water than in larvae or adults, irrespective of the origin of the sample. Moreover, larvae displayed a significantly different and most diversified microbial community than newly emerged adults, which appeared to be dominated by Proteobacteria. The microbiota of breeding site water significantly increased its diversity over time, suggesting the presence of a dynamic interaction among bacterial communities, breeding sites and mosquito hosts. The analysis of Wolbachia prevalence in adults from Foshan and five additional strains with different geographic origins confirmed the described pattern of dual wAlbA and wAlbB strain infection. However, differences in Wolbachia prevalence were detected, with one strain from La Reunion Island showing up to 18% uninfected individuals. These findings contribute in further understanding the dynamic interactions between the ecology of larval habitats and the structure of host microbiota, as well as providing additional information relative to the patterns of Wolbachia infection.

ViR: a tool to solve intrasample variability in the prediction of viral integration sites using whole genome sequencing data.

BACKGROUND: Several bioinformatics pipelines have been developed to detect sequences from viruses that integrate into the human genome because of the health relevance of these integrations, such as in the persistence of viral infection and/or in generating genotoxic effects, often progressing into cancer. Recent genomics and metagenomics analyses have shown that viruses also integrate into the genome of non-model organisms (i.e., arthropods, fish, plants, vertebrates). However, rarely studies of endogenous viral elements (EVEs) in non-model organisms have gone beyond their characterization from reference genome assemblies. In non-model organisms, we lack a thorough understanding of the widespread occurrence of EVEs and their biological relevance, apart from sporadic cases which nevertheless point to significant roles of EVEs in immunity and regulation of expression. The concomitance of repetitive DNA, duplications and/or assembly fragmentations in a genome sequence and intrasample variability in whole-genome sequencing (WGS) data could determine misalignments when mapping data to a genome assembly. This phenomenon hinders our ability to properly identify integration sites. RESULTS: To fill this gap, we developed ViR, a pipeline which solves the dispersion of reads due to intrasample variability in sequencing data from both single and pooled DNA samples thus ameliorating the detection of integration sites. We tested ViR to work with both in silico and real sequencing data from a non-model organism, the arboviral vector Aedes albopictus. Potential viral integrations predicted by ViR were molecularly validated supporting the accuracy of ViR results. CONCLUSION: ViR will open new venues to explore the biology of EVEs, especially in non-model organisms. Importantly, while we generated ViR with the identification of EVEs in mind, its application can be extended to detect any lateral transfer event providing an ad-hoc sequence to interrogate.

Population genomics in the arboviral vector Aedes aegypti reveals the genomic architecture and evolution of endogenous viral elements.

Horizontal gene transfer from viruses to eukaryotic cells is a pervasive phenomenon. Somatic viral integrations are linked to persistent viral infection whereas integrations into germline cells are maintained in host genomes by vertical transmission and may be co-opted for host functions. In the arboviral vector Aedes aegypti, an endogenous viral element from a nonretroviral RNA virus (nrEVE) was shown to produce PIWI-interacting RNAs (piRNAs) to limit infection with a cognate virus. Thus, nrEVEs may constitute a heritable, sequence-specific mechanism for antiviral immunity, analogous to piRNA-mediated silencing of transposable elements. Here, we combine population genomics and evolutionary approaches to analyse the genomic architecture of nrEVEs in A. aegypti. We conducted a genome-wide screen for adaptive nrEVEs and searched for novel population-specific nrEVEs in the genomes of 80 individual wild-caught mosquitoes from five geographical populations. We show a dynamic landscape of nrEVEs in mosquito genomes and identified five novel nrEVEs derived from two currently circulating viruses, providing evidence of the environmental-dependent modification of a piRNA cluster. Overall, our results show that virus endogenization events are complex with only a few nrEVEs contributing to adaptive evolution in A. aegypti.

Improved reference genome of the arboviral vector Aedes albopictus.

BACKGROUND: The Asian tiger mosquito Aedes albopictus is globally expanding and has become the main vector for human arboviruses in Europe. With limited antiviral drugs and vaccines available, vector control is the primary approach to prevent mosquito-borne diseases. A reliable and accurate DNA sequence of the Ae. albopictus genome is essential to develop new approaches that involve genetic manipulation of mosquitoes. RESULTS: We use long-read sequencing methods and modern scaffolding techniques (PacBio, 10X, and Hi-C) to produce AalbF2, a dramatically improved assembly of the Ae. albopictus genome. AalbF2 reveals widespread viral insertions, novel microRNAs and piRNA clusters, the sex-determining locus, and new immunity genes, and enables genome-wide studies of geographically diverse Ae. albopictus populations and analyses of the developmental and stage-dependent network of expression data. Additionally, we build the first physical map for this species with 75% of the assembled genome anchored to the chromosomes. CONCLUSION: The AalbF2 genome assembly represents the most up-to-date collective knowledge of the Ae. albopictus genome. These resources represent a foundation to improve understanding of the adaptation potential and the epidemiological relevance of this species and foster the development of innovative control measures.

Tracing temporal and geographic distribution of resistance to pyrethroids in the arboviral vector Aedes albopictus.

BACKGROUND: The arboviral vector Aedes albopictus became established on all continents except Antarctica in the past 50 years. A consequence of its rapid global invasion is the transmission of diseases previously confined to the tropics and subtropics occurring in temperate regions of the world, including the re-emergence of chikungunya and dengue in Europe. Application of pyrethroids is among the most widely-used interventions for vector control, especially in the presence of an arboviral outbreak. Studies are emerging that reveal phenotypic resistance and monitor mutations at the target site, the para sodium channel gene, primarily on a local scale. METHODS: A total of 512 Ae. albopictus mosquitoes from twelve geographic sites, including those from the native home range and invaded areas, were sampled between 2011 and 2018, and were analyzed at five codons of the para sodium channel gene with mutations predictive of resistance phenotype. Additionally, to test for the origin of unique kdr mutations in Mexico, we analyzed the genetic connectivity of southern Mexico mosquitoes with mosquitoes from home range, the Reunion Island, America and Europe. RESULTS: We detected mutations at all tested positions of the para sodium channel gene, with heterozygotes predominating and rare instance of double mutants. We observed an increase in the distribution and frequency of F1534C/L/S mutations in the ancestral China population and populations in the Mediterranean Greece, the appearance of the V1016G/I mutations as early as 2011 in Italy and mutations at position 410 and 989 in Mexico. The analyses of the distribution pattern of kdr alleles and haplotype network analyses showed evidence for multiple origins of all kdr mutations. CONCLUSIONS: Here we provide the most-up-to-date survey on the geographic and temporal distribution of pyrethroid-predictive mutations in Ae. albopictus by combining kdr genotyping on current and historical samples with published data. While we confirm low levels of pyrethroid resistance in most analyzed samples, we find increasing frequencies of F1534C/S and V1016G in China and Greece or Italy, respectively. The observed patterns of kdr allele distribution support the hypothesis that on site emergence of resistance has contributed more than spread of resistance through mosquito migration/invasions to the current widespread of kdr alleles, emphasizing the importance of local surveillance programs and resistance management.

Genomic Variant Analyses in Pyrethroid Resistant and Susceptible Malaria Vector, Anopheles sinensis.

Anopheles sinensis is a major malaria vector in Southeast Asia. Resistance to pyrethroid insecticides in this species has impeded malaria control in the region. Previous studies found that An. sinensis populations from Yunnan Province, China were highly resistant to deltamethrin and did not carry mutations in the voltage-gated sodium channel gene that cause knockdown resistance. In this study, we tested the hypothesis that other genomic variants are associated with the resistance phenotype. Using paired-end whole genome sequencing (DNA-seq), we generated 108 Gb of DNA sequence from deltamethrin -resistant and -susceptible mosquito pools with an average coverage of 83.3× depth. Using a stringent filtering method, we identified a total of 916,926 single nucleotide variants (SNVs), including 32,240 non-synonymous mutations. A total of 958 SNVs differed significantly in allele frequency between deltamethrin -resistant and -susceptible mosquitoes. Of these, 43 SNVs were present within 37 genes that code for immunity, detoxification, cuticular, and odorant proteins. A subset of 12 SNVs were randomly selected for genotyping of individual mosquitoes by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and showed consistent allele frequencies with the pooled DNA-seq derived allele frequencies. In addition, copy number variations (CNVs) were detected in 56 genes, including 33 that contained amplification alleles and 23 that contained deletion alleles in resistant mosquitoes compared to susceptible mosquitoes. The genomic variants described here provide a useful resource for future studies on the genetic mechanism of insecticide resistance in this important malaria vector species.

The Widespread Occurrence and Potential Biological Roles of Endogenous Viral Elements in Insect Genomes.

Modern genomic sequencing and bioinformatics approaches have detected numerous examples of DNA sequences derived from DNA and RNA virus genomes integrated into both vertebrate and insect genomes. Retroviruses encode RNA-dependent DNA polymerases (reverse transcriptases) and integrases that convert their RNA viral genomes into DNA proviruses and facilitate proviral DNA integration into the host genome. Surprisingly, DNA sequences derived from RNA viruses that do not encode these enzymes also occur in host genomes. Non-retroviral integrated RNA virus sequences (NIRVS) occur at relatively high frequency in the genomes of the arboviral vectors Aedes aegypti and Aedes albopictus, are not distributed randomly and possibly contribute to mosquito antiviral immunity, suggesting these mosquitoes could serve as a model system for unravelling the function of NIRVS. Here we address the following questions: What drives DNA synthesis from the genomes of non-retroviral RNA viruses? How does integration of virus cDNA into host DNA occur, and what is its biological function (if any)? We review current knowledge of viral integrations in insect genomes, hypothesize mechanisms of NIRVS formation and their potential impact on insect biology, particularly antiviral immunity, and suggest directions for future research.

Leaning Into the Bite: The piRNA Pathway as an Exemplar for the Genetic Engineering Need in Mosquitoes.

The piRNA pathway is a specialized small RNA interference that in mosquitoes is mechanistically distant from analogous biology in the Drosophila model. Current genetic engineering methods, such as targeted genome manipulation, have a high potential to tease out the functional complexity of this intricate molecular pathway. However, progress in utilizing these methods in arthropod vectors has been geared mostly toward the development of new vector control strategies rather than to study cellular functions. Herein we propose that genetic engineering methods will be essential to uncover the full functionality of PIWI/piRNA biology in mosquitoes and that extending the applications of genetic engineering on other aspects of mosquito biology will grant access to a much larger pool of knowledge in disease vectors that is just out of reach. We discuss motivations for and impediments to expanding the utility of genetic engineering to study the underlying biology and disease transmission and describe specific areas where efforts can be placed to achieve the full potential for genetic engineering in basic biology in mosquito vectors. Such efforts will generate a refreshed intellectual source of novel approaches to disease control and strong support for the effective use of approaches currently in development.

Polymorphism analyses and protein modelling inform on functional specialization of Piwi clade genes in the arboviral vector Aedes albopictus.

Current knowledge of the piRNA pathway is based mainly on studies on Drosophila melanogaster where three proteins of the Piwi subclade of the Argonaute family interact with PIWI-interacting RNAs to silence transposable elements in gonadal tissues. In mosquito species that transmit epidemic arboviruses such as dengue and chikungunya viruses, Piwi clade genes underwent expansion, are also expressed in the soma and cross-talk with proteins of recognized antiviral function cannot be excluded for some Piwi proteins. These observations underscore the importance of expanding our knowledge of the piRNA pathway beyond the model organism D. melanogaster. Here we focus on the emerging arboviral vector Aedes albopictus and we couple traditional approaches of expression and adaptive evolution analyses with most current computational predictions of protein structure to study evolutionary divergence among Piwi clade proteins. Superposition of protein homology models indicate possible high structure similarity among all Piwi proteins, with high levels of amino acid conservation in the inner regions devoted to RNA binding. On the contrary, solvent-exposed surfaces showed low conservation, with several sites under positive selection. Analysis of the expression profiles of Piwi transcripts during mosquito development and following infection with dengue serotype 1 or chikungunya viruses showed a concerted elicitation of all Piwi transcripts during viral dissemination of dengue viruses while maintenance of infection relied on expression of primarily Piwi5. Opposite, establishment of persistent infection by chikungunya virus is accompanied by increased expression of all Piwi genes, particularly Piwi4 and, again, Piwi5. Overall these results are consistent with functional specialization and a general antiviral role for Piwi5. Experimental evidences of sites under positive selection in Piwi1/3, Piwi4 and Piwi6, that have complex expression profiles, provide useful knowledge to design tailored functional experiments.